Soliton self-frequency blue-shift in gas-filled hollow-core photonic crystal fibers 
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We show theoretically that the photoionization process in a hollow-core photonic crystal hber hlled 
with a Raman-inactive noble gas leads to a constant acceleration of solitons in the time domain 
with a continuous shift to higher frequencies, limited only by ionization loss. This phenomenon 
is opposite to the well-known Raman self-frequency red-shift of solitons in solid-core glass fibers. 
We also predict the existence of unconventional long-range non-local soliton interactions leading 
to spectral and temporal soliton clustering. Furthermore, if the core is filled with a Raman-active 
molecular gas, spectral transformations between red-shifted, blue-shifted and stabilized solitons can 
take place in the same fiber. 

PACS numbers: 42.65.Tg, 42.65.-k, 32.80.Fb, 05.45.Yv, 42.81. Dp, 52.35.Sb 



Introduction — Hollow-core photonic crystal fibers (HC-PCFs) [l[ based on a kagome-lattice claddings have re- 
cently been shown to be very interesting for the investigation of broadband light-matter interactions between intense 
optical pulses and gaseous media. The fibers typically show transmission bands covering the visible and near-IR parts 
of the spectrum with relatively low loss and low group velocity dispersion (GVD), absence of surface modes, and high 
confinement of light in the core. Filled with a noble gas, they have recently been used in high-harmonic and efficient 
deep UV generation from fs pump pulses at 800 nm [H, Q . It has been previously shown that the Raman threshold 
can be drastically reduced in a HC-PCF filled with a Raman-active gas (such as H2) The system can be used for 
detailed experimental studies of, e.g., self-similar solutions of the sine-Gordon equation Q and backward stimulated 
Raman scattering 0, [1]. In bandgap-guiding gas-filled HC-PCFs, which have much narrower bands of transmission, 
a limited ionization-induced blue-shift of guided ultrashort pulses has been reported 0, Q . Very recently, ultrafast 
nonlinear dynamics in the ionization regime has been studied experimentally in Ar-fiUed kagome-style HC-PCF Q (a 
detailed account of these experiments is available in a parallel submission [10[). The reasons for the success of kagome 
HC-PCF in these applications are: (i) a group velocity dispersion (GVD) that is remarkably small (|/32| < 10 fs^/cm 
= 1 ps^/km from 400 to 1000 nm) in comparison to solid-core fibers (Fig. [TJa)) 0, [l^; (ii) the gas and waveguide 
contributions to the GVD can be balanced by varying the pressure, unlike in large-bore capillary-based systems where 
the normal dispersion of the gas dominates over the waveguide dispersion [llj . 

Photoionization in gases is traditionally modeled using the full electric field of the pulse [l2l |. In this paper, we first 
develop a new model to study pulse propagation in gas-filled HC-PCFs in terms of the complex envelope of the pulse. 
Using this model, we show analytically for the first time that intra-pulse photoionization leads to: (i) a soliton self- 
frequency blue-shift; (ii) long-range "non-local" soliton correlations and clustering; and (iii) spectral transformations 
of red-shifted, blue-shifted and stabilized solitons in Raman-active gas-filled HC-PCFs. 
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FIG. 1: (Color online), (a) GVD of an Ar-fiUed HC-PCF for gas pressures between 1 and 9 bar (calculated from [T3|). All 
subsequent calculations in this paper assume 5 bar pressure. Inset: cross-section of a broadband-guiding HC-PCF with a 
kagome-lattice cladding and a core diameter 30 /xm. Typical experimental transmission losses for the fundamental mode are 1 
dB/m for at 800 nm. (b) Comparison of the dependence of the Ar ionization rate on the pulse intensity using the full model 
of Eq. (m and the linearized model. 
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Governing equations — Photoionization can take place by either tunneling or multiphoton processes. These 



regimes are characterized by the Keldysh parameter pK l^il^. In the tunneling regime {pk ^ 1) the time- averaged 



ionization rate W(/) is given by [15|, |16[ 

= d {iH/lf' exp[-6 {lHllf'% (1) 

where d = 417o [ihY'"^ [Ui/UhV'^, b = 2/3 [Ui/Uh]^^'^, = 4.1 x 10^*^ Hz is the characteristic atomic frequency, 
Ui is the ionization energy of the gas 15.76 eV for argon). Uh ~ 13.6 cV is the ionization energy of hydrogen, 
Ih = 3.6 X 10^^ W/cm^ and / is the laser pulse intensity. For values of / in the range of 100 TW/cm^, the Keldysh 
parameter is pk ^ 1 for noble gases. However, experiments show that tunneling models provide excellent agreement 



with the experimental measurements even for ~ 1 [HI, ll8|- shown in Fig. Wlh), Eq- (Hlpredicts an ionization 



rate that is exponential-like for pulse intensities above a threshold value. Loss due to absorption of photons in the 
plasma is proportional to the ionization rate. Hence, any pulse with / 3> hh will have its intensity strongly driven 
back to near the threshold value, resulting in drastically reduced ionization loss. 

This allows us to use the first-order Taylor series to linearize the tunneling model just above / = Ith, where the 
optical pulses can survive for relatively long time without appreciable attenuation. Expanding Eq. ([T]) in its linear 
regime around an arbitrary point a = la/In, results in W w A/8 (A/), where A/ = / — /th, 5- = de~^ {2x — 
l)/[4a^/^///], Ith = alni'^x — b)/{2x — 1) is the threshold intensity, x — b/^/a, and a is chosen to reproduce the 
physically observed threshold intensity in the fiber of Fig. [ija), a = 2 x lO^'^. The purpose of the Hcavisidc function 
Q is to set the ionization rate to zero below the threshold intensity, see Fig. [Ijb). 

One can prove from first principles that propagation of light in a HC-PCF filled with an ionized Raman-active gas 
can be then described by the following coupled equations: 



id, + D{idt) + 7Ki?(t) ® - 7^ + *a 

ZKqC 



= 

, (2) 



where 5'(z,t) is the electric field envelope, z the longitudinal coordinate along the fiber, t is the time in a reference 
frame moving with the pulse group velocity, D(idt) = '^rn>2 l^rn{idt)™ / m\ is the full dispersion operator, (3^ is the 
TO-th order dispersion coefficient calculated at an arbitrary reference frequency wq, 7k is the Kerr nonlinear coefficient 
of the gas, R{t) = (1 — p)S{t) + ph{t) is the normalized Kerr and Raman response function of the gas, 5{t) is the 
Dirac delta function, p is the relative strength of the non-instantaneous Raman nonlinearity, h{t) is the causal Raman 
response function of the gas 0, the symbol ® denotes the time convolution, c is the speed of light, fcg = Uq/c, 
cjQ is the pulse central frequency, Wp = [e^nc/(eo"^c)]^^^ is the plasma frequency associated with an electron density 
nc(t), e and me are the electron charge and mass, and eq is the vacuum permittivity, a = ai -|- 0:2 is the total 
loss coefficient, ai is the fiber loss, a2 = dtn^ is the ionization-induced loss term, A^s is the effective mode 

area, Aj^fp = \^\'^ — l^'l^h, I^P = lA^a^ l^lth — -^th^oft, and nx is the total number density of ionizable atoms 
in the fiber, associated with the maximum plasma frequency wt = [e^nT/(eo'7ie)]^^^- In these coupled equations, 
the recombination process is neglected since the pulse duration (of the order of tens of fs) is always shorter than 
the recombination time [2^. If is measured in W, d/cAcs = 71 has the dimensions of W^^m^^. This is the 



nonlinearity associated with the plasma formation in the fiber. According to recent experimental measurements 2l|, 
7k shows a linear dependence on the gas pressure. These coupled equations ((2) are the first contribution of this paper. 
The validity of Eqs. © has been verified using a more complete ionization model based on the unidirectional wave 
equation [22| . 

Perturbation theory for floating pulses — In order to extract useful analytical information from Eqs. ([2|), further 
simplifications are necessary. For pulses with maximum intensities just above the ionization threshold (which we dub 
floating pulses, a new concept introduced in this paper for the first time), the ionization loss is not large and can 
be neglected as a first approximation. For such pulses, only a small portion of energy above the threshold intensity 
contributes to the creation of free electrons. Furthermore, for fioating pulses one can remove the 8-function from 
the equations, provided that the cross-section a is replaced by a properly reduced a' that takes into account the 
overestimation of the ionization rate (23j . Introducing the following rescalings and redefinitions: ^ = z/zq, t = t/to, 
*o = [iKZo]^^^'^, tp = */^'o, r{T) = R{t)to, (j) = ifcozo [^p/^o]^, 0t = ^^o^o [wt/wq]^, and a = ct' io/[^off7K ^o], 
where Zq = tQ/|/32(wo)| is the second-order dispersion length at the reference frequency cjq and t^ is the input pulse 
duration [Toj . Hence, the two coupled equations for floating pulses can be replaced by 



id^ + D{idr)+r{T)(S>\ip{T)\^ 



(3) 
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The total number of photons is conserved in this set of coupled equations - in contrast to Eqs. ([2]) - since losses are 
neglected for floating pulses. 

The effect of the Raman and ionization perturbations on the soliton dynamics in HC-PCFs can be studied using 

Eqs. ([3]). The second equation can be solved analytically, 0(r) = 0t |l — cxp —a \'tp{T')\'^dT' |, with the initial 

condition 0(— oo) = 0, corresponding to the absence of any plasma before the pulse arrives. For a small ionization cross- 
section, 0(r) ~ 77 |?/'(T')pdr', where 77 = ct^t- Moreover, in the long-pulse limit |V'(t— t')P ~ \ip{T)\'^ — T'dT\ilj{T)\ 
[l9| . This allows the two coupled equations to be reduced to a single partial integro-differential equation: 

id^ij + D{idr)ij + Ki V - TR.V'^rlV'P -Vi^ f \tp\^dT' = (4) 

<y —00 

where tr = t' r{T') dr' . This equation shows clearly that the effect of ionization is exactly opposite to that 
of the Raman effect: the fourth term in Eq. ^ involves a derivative of the field intensity, while the fifth term 
involves an integral on the same quantity. One can then conjecture that the last term will lead to a soliton self- 
frequency blue-shift due to ionization, instead of a red-shift. To prove this statement, we use the perturbation 
theory described in [l^. First, the soliton functional shape is assumed to be unchanged during the action of the 
perturbations induced by the Raman effect and the photoionization process (this must be verified a posteriori): 
''Psi£,,T) = Aosech[AQ{T ~ Tp{S,))] e~'^^^^^'^ , with Tp(^) is the temporal location of the soliton peak and is the 

self-frequency shift. When this Ansatz is inserted into Eq. (|4]), simple ordinary differential equations can be obtained 
for both n{^) and rp(^), results in = f2Raman(6 + ^^ion(C) = -9^, Tp{£,) = jl, and g = grcd + ffbiuo, where 
5icd = +(8/15)rF.AQ and gbiuo = — (2/3)?7^q. Note that g can be positive, negative or even zero, depending on 
the value of -q, tr and Aq. By using the exact solution for (/)(r) given previously, one obtains the more precise 
rate g{^i^g = cr~^Ag ^(^t [(1 — cfAq) — {1 + ctAq) exp(— 2crAo)], which tends to gune for small values of a, but starts 
to differ considerably from it for Aq > a^^ . The above solution clearly shows that, in the range of validity of 
perturbation theory (i.e., for floating solitons), "photoionization leads to a soliton self-frequency blue-shift. This blue- 
shift is accompanied by a constant acceleration of the pulse in the time domain - opposite to the Raman effect, which 
produces pulse deceleration. The blue-shift is limited only by ionization loss, which slowly decreases the pulse intensity 
until it falls below the threshold value. Other important effects, such as the formation of inverted gravity-like bound 
solitons supported by the plasma nonlinearity, will be reported elsewhere. 

In the presence of ionization-induced losses above the threshold intensity, Eqs. ([2]) must be numerically solved to 
study the full dynamics of floating pulses. Figures [Ha,b) show the temporal and spectral evolution of a high-order 
input soliton, closely following the results reported in the companion experimental paper [lo| . When the intensity of 
the energetic pulse exceeds the threshold value as a result of self-compression, a fundamental soliton is ejected from 
the main pulse and continues to blue-shift until ionization loss reduces its amplitude below the threshold value. At 
longer distances, another compression occurs and a second soliton is generated. The use of a kagome-style HC-PCF is 
essential to observe the soliton blue-shift, since conventional photonic-bandgap flbers have much stronger dispersion 
variations, which would quickly destabilize any possible solitary wave as in [7|, |8|. 

Long-range non-local soliton forces and clustering — An interesting and unexpected interaction occurs between 
two solitons when their temporal separation is shorter than the recombination time, due to the non- vanishing electron 
density tail. Using the exact formula for the ionization field </>(t), one can see that a leading soliton with amplitude 
Aq can slow down the acceleration of a trailing soliton by an exponential factor exp(— 2crAo). The reason is that the 
ionization field 4>{t), created by the first soliton, decays at a relatively slow rate. This establishes a unique non-local 
interaction between this soliton and other temporally distant solitons. 

Figures [Ua,b) show the output temporal and spectral dependence of a pulse iVsechr on the soliton order TV. In 
the presence of ionization loss, when the intensity of the leading soliton decreases to the threshold value, i.e., the 
blue-shifting process ceases, the lagging soliton will recover its expected blue-shift. The reason is the disappearance 
of the exponential decaying factor at that particular point. Also, there is a maximum frequency attained by each 
soliton that depends on the initial soliton intensity. These interactions may lead, at some 'magic' input energy, to 
clustering two or more distinct solitons in both temporal and spectral domains, as shown in Figs. [3l^a,b). 

Soliton spectral transformations — Interestingly, this perturbation theory for floating solitons predicts the forma- 
tion of spectrally stabilized solitons in Raman-active gases due to the different signs and ^o-dependence of the Raman 
and photoionization shifts, g^ed oc Aq (a well-known result [Toj). and gbiue oc —Aq (reported for the first time in this 
paper). If one launches an sufficiently energetic pulse into the fiber, soliton fission takes place, independent of the 



particular perturbation applied [24|, l25| . This generates a train of fundamental solitons with progressively decreasing 
peak amplitudes. 

The temporal and spectral evolution of such a pulse when it propagates in a mixture of argon and air (Raman-active) 
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FIG. 2: (Color online). Temporal (a) and spectral (b) evolution of an energetic pulse propagating in the Ar-fiUed HC-PCF. 
The temporal profile of the input pulse is A^sechr, with N — 8. The panels show the ejection of two fundamental solitons that 
continuously blue-shift until ionization loss reduces their intensities below the threshold value. 
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FIG. 3: (Color online). Temporal (a) and spectral (b) outputs of an energetic pulse A'^sechr after propagating inside an 
Ar-filled HC-PCF with length ^ = 1/4 versus the soliton order A^. Temporal and spectral clustering occur at N = 9.2 due to 
the long-range "non-local" soliton interactions described in the text. 
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is depicted in Figs. |4l^a,b). After the fission process, solitons with intensities less than the threshold intensity are 
red-shifted by the undisturbed Raman process. However, solitons possessing intensities above the threshold value are 
influenced simultaneously by both the photoionization and the Raman effects. Depending on their initial intensities, 
these solitons can be initially blue-shifted, red-shifted or stabilized. The Raman self-frequency red-shift will be more 
pronounced than the ionization self-frequency blue-shift for floating solitons possessing initially larger amplitudes 
{Aj > Acr), where Aa- is a critical amplitude. However, for less intense floating solitons {Aj < Aa) it may happen 
that exactly the opposite phenomenon occurs, i.e., the blueshift will dominate. The critical intensity can be estimated 
from the equation g^cd + 5biuc = 0, giving A^j. = 5r]/(ATYi). When the ionization loss arrests the photoionization 
process, the initially blue-shifted solitons start to reverse their self-frequency shift towards the red. At a certain point, 
these solitons can become frequency-stabilized over a short distance. This may result in multiple collisions between 
floating solitons if their temporal trajectories intersect. When the instantaneous intensity exceeds the threshold value 
upon collision, a second blue-shift event may occur. 

Conclusions — A direct photoionization process can act on solitons by constantly blue-shifting their central 
frequencies representing the exact counterpart of the Raman self-frequency red-shift when the intensity of solitons 
is slightly above the photo-ionization threshold. This spectral transformation is limited by the ionization loss that 
restricts the pulse intensity to the threshold value, hence arresting the soliton blue-shift. The new theoretical model, 
presented by Eqs. ([2]), is suitable for analytical manipulations, and has led us to predict a number of new phenomena 
such as long-range non-local correlation forces, and spectral transformation between red- and blue-shift in Raman- 
active gases. The results reveal new physics and offer novel opportunities for the manipulation and control of the 
soliton dynamics inside these versatile optical waveguides. 
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FIG. 4: (Color online). Temporal (a) and spectral (b) evolution of an energetic pulse propagates in a HC-PCF filled with argon 
and air. The temporal profile of the input pulse is A''sechr, with A'' = 4. Soliton temporal and spectral trajectories show an 
initial acceleration and blue-shift due to plasma formation, an intermediate stabilization below-threshold intensity, and finally a 
deceleration and red-shift due to Raman effect. Near the fiber end a second blue-shift event takes place due to soliton collision, 
generating a second surge of ionized plasma. 



